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Abstract 

We present results obtained for nuclei, nuclear and neutron star matter, 
' and neutron star structure obtained with the recent Argonne vis two- 

\^ ■ nucleon and Urbana IX three-nucleon interactions including relativistic 

boost corrections. These interactions predict that matter will undergo 
a transition to a spin layered phase with neutral pion condensation. We 
■ also consider the possibility of a transition to quark matter. 

00 ' 

O ■ 1 Introduction 

<N , 

Properties of matter having densities up to ~ 1 nucleon/fm 3 are important 
in determining the structure and maximum mass of neutron stars. At such 
densities the average interparticle spacing is ~ 1 fm. It is well known that 
the rms charge radius of proton is ~ 0.8 fm, thus one may be concerned that 
nucleons 1 fm apart have too large a structural overlap, and may even cease 
to be nucleons. However, the large charge radius of protons is due to a diffuse 
cloud attributed to mesons. The charge form factor of the proton, as well as 
the magnetic form factors of the proton and the neutron are well approximated 
rS [ by the dipole form: 



=5 



F(q) = (l + q 2 /q 2 )- 2 ; q = 0.84 (1) 

The charge density obtained by inverting this form factor is: 

p ch {r) = 3.3 e~ r/0 - 23 fm" 3 , r in fm. (2) 

The charge densities of two protons placed at ±0.5 fm on the z-axis are shown 
in fig. 1. We see that there is not too much structural overlap of protons 1 
fm apart. We may therefore assume that they still are protons, and absorb 
the effects of the overlap into the two nucleon interaction. This assumption is 
supported by studies of the deuteron. The deuteron wave function peaks at 
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Figure 1: Charge density of two nucleons at z — ±0.5 fm, along the z-axis 

~ 1 fm separation between the neutron and the proton in agreement with the 
observed deuteron form factors 

In nuclear many-body theory, nuclei and nuclear matter are described by 
the Hamiltonian: 



in which i,j, k, ... denote nucleons. The first term is the nonrelativistic kinetic 
energy, and the second contains the two nucleon interaction Vij in the center of 
mass frame of the interacting pair. The relativistic boost interaction 5v(Pij) 
describes the dependence of two-nucleon interaction on their total center of 
mass momentum Py, and the last term contains the three nucleon interaction. 
In principle the H can have additional terms, such as the boost correction for 
three-nucleon interaction, four nucleon interaction, etc.; however these terms 
are believed to have a negligible effect on neutron star structure. The above 
H describes only the hadronic part of neutron star matter which also contains 
electrons and, at high densities, muons to preserve charge neutrality and beta 
equilibrium. The relativistic kinetic energies of the leptons and the Coulomb 
interaction energies are to be added to the above hadronic H. 

Recent developments in various terms of the nuclear Hamiltonian are dis- 
cussed in the next section, while progress in the variational calculations of the 
properties of matter from this H is reviewed in sect. 3. The results for neutron 
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stars are given in sect. 4, while the possibility of a phase transition to quark 
matter is considered in sect. 5. A more detailed description of this work will 
be published elsewhere. 

2 Modern Models of Nuclear Forces 

In the early 1990's the Nijmegen group || carefully examined all the data 
on NN scattering at energies below 350 MeV published between 1955 and 
1992. They extracted 1787 proton-proton and 2514 proton-neutron "reliable" 
data, and showed that these could determine all NN scattering phase shifts 
and mixing parameters quite accurately. NN interaction models which fit this 
Nijmegen data base with a x 2 /^data ~ 1 ar e termed "modern". These include 
the Nijmegen models || called Nijmegen I, II and Reid 93, the Argonne Vi& Q 
(A18) and CD-Bonn ||. In order to fit both the proton-proton and neutron- 
proton scattering data simultaneously and accurately, these models include a 
detailed description of the electromagnetic interactions and terms that violate 
the isospin symmetry of the strong interaction via the differences in the masses 
of the charged and neutral pions, etc. 
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Figure 2: The E(p) of neutron matter calculated from modern potential mod- 
els. The upper four curves show results obtained with Brueckner-Hartree-Fock 
calculations, while the lower curve is with variational calculations. 

In these five models the two nucleon interaction is expressed as a sum 
of the one-pion exchange potential vfj and the rest of the interaction vfy They 
use different parameterizations of the v^, and the Nijmegen-I and CD-Bonn 




also include nonlocalities suggested by boson-exchange representations. Thus, 
like the older models, they make different predictions for many-body systems. 
However, the differences in their predictions are much smaller than those be- 
tween older models, presumably because they exactly fit the same large data 
set. For example, the triton ground state energies predicted by the modern 
Nijmegen and Argonne models, without the 5v(Pij) and Vijk terms in the 
Hamiltonian, are between -7.62 and -7.72 MeV ||, while that of CD-Bonn 
is -8.00 MeV |J. This should be compared with the spread from -7.35 MeV 
(original Reid) to -8.35 MeV (Bonn-A) in the predictions of the older models. 
The predictions for the energies of dense neutron matter obtained from the 
modern potentials using lowest order Brueckner-Hartree-Fock (LOB) method 
[0] are also quite close together as shown in fig. 2. This figure also shows the 
energies of neutron matter calculated from the A18 potential by the variational 
method using chain summation techniques j5|. The difference between the re- 
sults obtained for this modern potential using LOB and variational methods 
is larger than that in the results of all the modern potentials with the LOB 
method. 

The modern potentials, like their predecessors, do not explain all of the 
observed nuclear properties. They all underbind the triton, and overpredict 
the density of nuclear matter. It is necessary to add three-nucleon (NNN) 
interactions to the Hamiltonian to reproduce these observables. The Urbana 
models of NNN interaction have only two terms: 

V ijk = V% + (4) 

the first gives the two-pion exchange NNN interaction via the pion-nucleon 
delta resonance 0, and its strength is denoted by the parameter A 27T ; the sec- 
ond is a phenomenological spin-isospin independent, shorter range term with 
strength denoted by Uq. The values of the parameters A 2n and Uo of the 
Urbana-IX (UIX) NNN potential are chosen such that in combination with 
A18 NN potential triton energy is reproduced via exact calculations and the 
equilibrium density of nuclear matter is reproduced via variational calcula- 
tions PH[| . The boost interaction 8v(Pij) was negelected in these calculations 
to determine A 2n and Uq. Its effect is discussed later. It is unlikely that a 
combination of CD-Bonn and UIX potentials will reproduce the triton en- 
ergy. Parameters of the NNN interaction have to refitted for each model of 
the NN interaction, as in the recent work of the Bochum group [Tl| using the 
Tucson- Melbourne |12| model of VM. The fitted Vijk then partly corrects for 



the deviation of the model Vij from its exact representation. We hope that 
predictions of the modern combinations of and V^k will be significantly less 



Table 1: Results of Quantum Monte Carlo Calculations in MeV 
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model dependent than those of the i>jj alone. 

It has recently become possible to calculate all the bound states of up to 
eight nucleons from realistic nuclear forces with the Greens Function Monte 
Carlo (GFMC) method. Since solar and primodial fusion reactions primarily 
involve nuclei having A < 8, they have a special role in the universe. Here we 
use the results obtained with A18, without 5u(Py), and UIX interactions to 
test the accuracy of variational wave functions and the interaction models. The 
calculated energies from [ID] and [13 are listed in Table I. Its first two columns 
specify the nuclear state and the next four give the calculated expectation 
values of the interaction components. The last three columns list the energy 
calculated with GFMC, the difference between the experimental and GFMC 
energies, and that between our optimum variational and GFMC energies. The 
v ^ seems to give the dominant contribution to nuclear binding, while v R is also 
essential for all but the deuteron. There is a large cancellation between the 
kinetic and two-nucleon interaction contributions, causing the total nuclear 
energy to be much smaller than the (vij). The (V^£) is much smaller than 
(v^); however, it is a significant fraction of the nuclear binding energy. The 
(V^.) is the smallest, and AE expt ., the difference between experiment and 
theory is even smaller. In the neutron rich nuclei 7 He and 8 He, the AE expt . is 
comparable to the (V^), indicating that the UIX model may not be describing 
the interaction between three neutrons very well. All the p-shell nuclei having 
A > 5 are underbound, and the problem increases with the magnitude of 
nuclear isospin T. New models of are being studied to reduce the AE expt ,. 

Table 1 also reveals another problem; the difference AEvmc between GFMC 



and variational Monte Carlo (VMC) energies is surprisingly large in the p-shell 
nuclei. Earlier successes in variational calculations of the A < 4 s-shell nuclei 
led to hopes that the variational calculations of larger nuclei and uniform nu- 



cleon matter, using cluster expansions |H |, and chain summation methods || 
may have less than 10 % errors. The present VMC for 8 Be, which includes 
many three-body correlations omitted in variational calculations of larger nu- 
clei and uniform matter, has a 12 % error. It appears that important aspects 
of the wavefunctions of p-shell nuclei are still not understood. 

All NN interaction models are obtained by fitting NN scattering data in 
the center of mass frame. The model Vij denotes the interaction between two 
nucleons in the frame in which their total momentum Py = pj + pj, is zero. 
In the rest frame of all nuclei other than the deuteron, the Py of a pair of 
nucleons is non-zero. One then has to use the correct NN interaction: 

v(P ij ) = v ij + 5v(P ij ), (5) 
between particles with total momentum P#. The correction 5v(Pij) is called 



the boost interaction [p~5[]; it is zero when P^- = 0. 



It is useful to consider a familiar example. The Breit interaction |L6] be- 
tween two particles of mass m and charge Q is, ignoring spin dependent terms 
for brevity, given by: 




Pi ' Pj Pi ' *ijPj ' 



2m 2 2m 2 r. 2 ; 



(6) 



which depends upon both pj and pj. We can express it in our notation as a 
sum of and 5v(Pij): 

13 r tJ { + 2m 2 + 2m 2 r% J' UJ 

^ = - % G§ + ^Bf) • (8 > 

where pjj = (p« — Pj)/2, is the relative momentum. The dependence of on 
Pij is included in all modern models, however, the 5v(Pij) has been neglected 
in the majority of nuclear and neutron star calculations. 

Following the work of Krajcik and Foldy |L7], Friar [18| obtained the fol- 



lowing equation relating the boost interaction of order P 2 to the interaction 
in the center of mass frame: 

MP) = + A[P • rP ■ V, v] + -Lfto - a 3 ) x P • V, i;]. (9) 



The general validity of this equation in relativistic mechanics and field the- 
ory was recently discussed Hl5| . Nonrelativistic Hamiltonians containing boost 
interactions include all terms quadratic in the particle velocities. The con- 
tribution of the two-body boost interaction to the energy of light nuclei has 
been evaluated with VMC method fll^, |2(], 21 1; it is repulsive and equals ~ 37 
% of that of V# k listed in Table 1. Therefore about 37 % of the Vg k in UIX 
simulates the contribution of the neglected Sv. The three- nucleon interaction 
to be used in Hamiltonians containing Sv is denoted by V£ k ; the strength of 
V*$ is 0.63 times that of V { f k in UIX, while V$? = Vgj. In nuclei and in 
nuclear matter at densities up to ~ 0.2 nucleons/fm 3 Hamiltonians containing 
V* + Sv and V alone give rather similar results, but at higher densities they 
differ substantially. Naturally the results of the Hamiltonian with V* + Sv are 
more reliable. 

One can also consider relativistic nuclear Hamiltonians of the type: 



(10) 



where v, V, ... include relativistic nonlocalities 21 , and the boost interactions 



include terms of higher order in Pj,-. We must refit the two-nucleon scatter- 



19, 21 



ing data to determine the v in Hr, using relativistic kinetic energies 
Uncorrelated nucleons in nuclei have small momenta of order m/4, thus the 
correction to their nonrelativistic kinetic energy is negligible. However, due 
to correlations induced by the strong v, a pair of nucleons acquire large rela- 
tive momenta at small r^-. The results obtained for light nuclei with Hr and 



the H given by eq.(|3) are very similar pi|, indicating that substantial im- 



provements are obtained by including the relativistic boost interactions in the 
nonrelativistic Hamiltonian, as is well known for electromagnetic interactions. 



3 Dense Nucleon Matter 

The equation of state (EOS) of cold symmetric nuclear matter (SNM) and 
pure neutron matter (PNM) has been recently calculated with the variational 



method using chain summation techniques || |22|. Calculations have been 
done with the Hamiltonians with and without boost interactions. 

We will first discuss the results presented in ref. || for the Hamiltonian 
without boost interactions. The calculated SNM and PNM energies are shown 



in figs. 3 and 4 along with the results of earlier calculations [^3|, 2j]. The density 
dependence of the U-DDI interaction was choosen to obtain the empirical 
saturation properties of SNM, so the minimum of that curve may be regarded 
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Figure 3: E(p) of SNM calculated from A18 and UIX interactions compared 
with the results of earlier calculations. The two sets of variational minima 
obtained at p > 0.28/ fm 3 are labeled LDP and HDP for low and high density 
phases. 



as experimental data. The calculated energy of —12 MeV/A is higher than the 
observed —16 MeV/A. However we now believe that much of this difference is 
due to the simplicity of the variational wave functions used in nucleon matter 
calculations. Use of improved variational wave functions should lower the SNM 
energy by more that two MeV/A. Comparison of VMC and GFMC results for 
energies of eight neutrons bound in a weak potential well |25| indicate that the 
present variational wave functions are more accurate for PNM than for SNM. 

The calculated energies indicate a phase transition in SNM at a density 
of ~ 0.3/ fm 3 (fig.3) and in PNM at ~ 0.2/ fm 3 (fig.4). Detailed analysis of 
the pair distribution functions and sums of response functions || indicate that 
the phase at higher densities has spin-isospin order expected from neutral pion 



condensation, first considered by Migdal ||26|| . The variational wave functions 



used in this study are appropriate for isotropic matter. The phase transition 
is signaled by a change in the range and strength of tensor correlations. It is 
likely that the phase at higher density will be better described with a liquid 
crystal wave function with alternating spin layers P7| used in many studies 
with effective interaction models. 

As realized long ago by Migdal, this transition is very sensitive to the short 
range behavior of the <7j • <jj and cr, • OjTi ■ Tj interactions between nucleons. 
It did not occur in either SNM or PNM with the Urbana Vu model of 1981, 
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Figure 4: E{p) of PNM. See fig. 3 caption for details. 

while it occured only in PNM with the later Argonne v\± model of 1984. The 
Urbana-Argonne potentials have similar forms, but are fit to different data sets. 
The 1981 (1984) Urbana (Argonne) vn models were fit to the n-p phase shifts 
available in the late 1970's (early 1980's), while the A18 is fit directly to the 
Nijmegen 1994 p-p and n-p scattering data base. The quality of the fit obtained 
by the A18 potential is much higher, and thus it is likely that it provides a 
more accurate representation of the nuclear force. The pion exchange part, 
V?l, of the three- nucleon interaction is essential for the transition to occur in 
SNM, while in PNM it occurs at a much higher density in the absence of V^. 

The composition of nucleon matter energy calculated from A18 and UIX 
interactions is listed at selected densities in table 2, in which T-1B denotes the 
Fermi gas kinetic energy. The V-2B-S and T-2B-S list the two-body (2B) clus- 
ter contributions to the potential and kinetic energy from the static parts of the 
pair interaction and correlation operators. The V-2B-MD and T-2B-MD show 
2B contributions having one or more momentum dependent (MD) interaction 
or correlations. These are conveniently regarded as the difference between the 
total 2B energies and their static parts. The static many-body (MB) contribu- 
tions, MB-S, as well as the expectation values of the three- nucleon interactions 
are calculated with the Fermi hypernetted and single operator chain summa- 
tion methods, presumably quite accurately. However, the MB contributions 
containing MD interaction or correlations, as well as the boost interaction 
5v(Pij), are evaluated from dressed three-body clusters. The 5E-2B lists an 
estimate of the change in energy that may be obtained by allowing additional 




Table 2: Composition of nucleoli matter energy in MeV. 
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flexibility in the two-body correlations. We note that at high densities the MD 
parts, the V*j k and 5v(Pij) give significant contributions. 

4 Neutron Star Matter 

The matter at subnuclear densities in the outer and inner crusts of neutron 
stars has interesting structures as discussed at this meeting by Pethick and 



Ravenhall |28|, |29j. Here we focus on the matter below the crust assuming 
that it is a uniform mixture of neutrons, protons, electrons and muons in 
beta equilibrium. The calculated energies of SNM and PNM are fitted by 



generalized Skyrme type effective interactions |3(| having different parameters 
below and above the pion condensation phase transition. The beta equilibrium 
conditions are calculated from these effective interactions. 

Ignoring mixed phase regions the normal matter at density and proton frac- 
tion of 0.204 fm~ 3 and 0.073 is found to be in equilibrium with pion condensed 
matter at 0.237 fm~ 3 and 0.057 for our most reliable model with boost and 
three-nucleon interactions. Obviously the matter in the phase with pion con- 
densation has a lower charge density than the matter without condensation. 
Therefore, in reality the transition will proceed through mixed phase regions 



of the type discussed by Glendenning |T| and Heiselberg et. al. |32| in the 



context of the transition from hadronic to quark matter. The mixed phase re- 
gions in the pion condensation transition do not seem to have a large effect on 



the structure of neutron stars because the discontinuity in the charge density is 
rather small. For example, the predicted thickness of the mixed phase regions 
varies from ~ 40 to 14 m in stars with 1.41 to 2.1 M ; however these predic- 
tions are quite crude because our calculations of the pion condensed phase are 
still incomplete. 

The calculated density dependence of the proton fraction in neutron star 
matter is shown in fig. 5 for the Hamiltonians containing only A18, A18+ 5v , 
Al8+5t> +UIX*, and A18+UIX interactions. The earlier results obtained with 
the U-DDI interaction (FPS) are shown for comparison. Both the boost and 
the three nucleon interactions increase the proton fraction in matter, but it 
remains below the critical value of 0.148 needed for direct Urea cooling, in 
the range of neutron star densities. The plus signs in fig. 5 show the proton 
fractions calculated with the A18 interaction alone with the LOB method 0. 
These are in reasonable agreement with the results of our calculations up to a 
density of 0.6 fm -3 . The density dependence of the electron chemical potential 
in matter is shown in fig. 6. 
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Figure 5: The proton fraction in matter in beta equilibrium for various model 
Hamiltonians. The plus signs show the results obtained for the A18 Hamilto- 
nian with LOB method. 

The predicted neutron star properties appear in figures 7 and 8. The 
dotted lines in these figures show results obtained using the PNM EOS. They 
are not very different from those obtained with matter in beta equilibrium. 
Due to the presence of momentum dependent and three-nucleon interactions 
in the present nuclear Hamiltonians, the predicted sound velocity in matter 
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Figure 6: The electron chemical potential in neutron star matter for various 
model Hamiltonians. Matter contains muons when /z e exceeds the muon mass 



can exceed the velocity of light. The densities at which super luminal sound 
occurs in each model are marked by vertical bars in fig. 7. These densities are 
very close to the maximum densities that can occur in neutron stars, therefore 
limiting the the velocity of sound to be < c does not have a significant effect 
on the predicted maximum masses. 

Observations of binary neutron stars have confirmed the existence of 1.4 
M & neutron stars allowed by all models. Recently several authors [| 



34, 35 



have argued that there are indications of the existence of neutron stars with 
M ~ 2M , which would rule out models without three-body forces, but see 
ref. |36| and Lamb's contribution p7|. 



5 Transition to Quark Matter 

It is expected that at some large density there will be a transition from nucleon 
matter to quark matter (QM). Several authors jnj have studied the energy of 
cold QM using the Bag-model, in which the total energy density contains a 
"Bag-constant" B, that takes into account the difference between the energies 
of the physical and QM vacua, and the energy of quarks interacting via one 
gluon exchange interaction calculated in first order of a s . The u and d quarks 
are assumed to be massless, and mass of s quarks is taken as 150 MeV. The 
energy density of QM obtained with B = 122 and 200 MeV/fm 3 is compared 




Figure 7: The dependence of neutron star mass on its central density for various 
model Hamiltonians. The full (dotted) lines show results obtained with the 
EOS of matter in beta equilibrium (PNM), and the vertical bars show where 
matter becomes superluminal in these models. 




Figure 8: Neutron star mass-radius relation obtained from various model 
Hamiltonians. See caption of fig. 7 for notation. 



with that of nucleoli matter (NM) in beta equilibrium for the more realistic 
models containing boost interaction, in fig. 9. The value B = 122 MeV/fm 3 
is supported by an analysis using the average of nucleon and delta-resonance 
masses |39| |. In the absence of three-body interactions NM is found to have 
lower energy than QM up to a large density of ~ 1.6 fm -3 . However, the 
transition density is lowered to ~ 1 fm -3 after including the contributions of 
the three nucleon interaction. 
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Figure 9: The energy densites of electrically neutral QM and NM are shown 
by full lines, while the dashed lines show those for matter with mixed QM and 
NM phases. 



Glendenning |3T| realized that it is not necessary to require the QM and 
NM phases to be separately charge neutral; a mixture of NM and QM can exist 
in a uniform lepton gas. Neglecting the energy of the interface between the 
QM and NM, and using the A18 + 5v + UIX* model of NM, the transition 
occurs over the density range p = 0.74 to 1.80 fm -3 for B = 200, and 0.58 
to 1.46 fm~ 3 for B = 122 MeV/fm 3 . At the lower end of this density range 
we have mostly NM with drops of QM, while at the higher end we will have 
mostly QM with drops of NM p2|. 



The dependence of the mass of neutron stars on their central density, ob- 
tained after incuding the effects of the transition to quark matter on the EOS 
are shown in fig.10. For B = 122 Mev/fm 3 and the A18 + Sv + UIX* NM 
Hamiltonian stars with masses above 1.5 M seem to have drops of quark mat- 
ter in their cores. The maximum mass is reduced to ~ 1.9M Q corresponding 
to a central density of ~ 1.1 fm~ 3 . In this model pure quark matter appears at 
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Figure 10: Dependence of neutron star mass on its central density. Full lines 
show results obtained with NM EOS, while dashed lines show results obtained 
after including the effect of the transition to QM on the EOS. 

p = 1.46 fm -3 ; therefore even the most massive stars have mixed QM and NM 
phases in their interior. On the other hand, if B = 200 Mev/fm 3 the maximum 
mass becomes ~ 2M Q , and stars with masses from 1.93 to 2.01 M have quark 
drops in their interior. 
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